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Abstract 

We evaluate the contribution of second order terms to the pion- 
nucleus s-wave optical potential of pionic atoms generated by short 
range nuclear correlations. The corrections are sizeable because they 
involve the isoscalar s-wave vrA^ amplitude for half off-shell situations 
where the amplitude is considerably larger than the on-shell one. 

In addition, the s-wave optical potential is reanalyzed by looking at 
all the different conventional contributions together: lowest order, Pauli 
corrected rescattering term, second order absorptive effects, terms from 
the interaction of pions with the virtual pion cloud (chiral corrections) 
and correlation effects. Different off-shell extrapolations for the ttN 
amplitude are used and it is found that, although some individual terms 
are sensitive to the extrapolation, the sum of them is rather insensitive. 
The results are compared with empirical values from best fits to the 
data and are found to be compatible, within theoretical and empirical 
uncertainties. The results do not rule out further contributions but they 
put stringent constraints on their strength. 
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The original work of Ericson-Ericson on pionic atoms [Q set already the 
two largest contributions to the s-wave pion nucleus optical potential. Its real 
part is given by 

RenW = -47r[(l + e)(6o + A6o)p + (1 + e)6i(p„ - Pp) + (1 + ^e)Re5op'] (1) 
with e = f^, p = p„ + Pp and 

&o = ^(fli + 203) ; 61 = -^(01-03) (2) 

where Oi, 03 are the isospin |, | scattering lengths given by 

ai = (0.171 ±0.004) m;^; 03 = (-0.105 ± 0.003) m^^ (3) 

The magnitude A60, introduced in [jl| and subsequently rederived in differ- 
ent formalisms with exactly the same results ^, accounts for the Pauli 
corrected second order rescattering terms. Its derivation within the many- 
body framework used in the present paper is shown in appendix A of ref. |Q 
and its value is given by 

ah. 1 1 

Ah, = -^-—{\l + Xl)- Ai,2 = --m.(l + e)6o,i (4) 
■nml 1 + e 2 

In addition, there are other terms of higher order in the nuclear density, 
accounted for in terms of ReSo, which in ref. [IH were fit to the data, but 
which have been subsequently evaluated. We discuss them below. 

One of them is the real part corresponding to the absorption diagrams, 
also known as dispersive correction, which is calculated in [|, |[. The results 
were obtained using two different off-shell extrapolations, one from Hamilton 
and the other one from LMM [^f. The uncertainties from different sources 
were estimated to be of the order of 30%, the largest one coming from the off- 
shell extrapolation. Indeed, the results obtained there in terms of the second 
order parameter Re-Bo were 

5Re5o,d = 0.032 m'^ [5] ; 5Re5o,rf = 0.017 m"^ [6] (5) 

using the off-shell extrapolations of ref. |]^ (Hamilton) and ref. |^ (LMM) 
respectively. 

The findings in ref. on the effective density peg = 0.5 po (po = —0.17 fm~^) 
for the s-wave part of the optical potential, allow one to talk alternatively in 
terms of an equivalent 56o,d parameter, multiplying (5-Bo,d by 0.23 m^. 

On the other hand, recent evaluations of the contribution to the s-wave 



self-energy from the interaction of pions with the virtual pion cloud [|Ty, |TT 



within the context of Weinberg chiral Lagrangians |T^, gave rise to a moderate 
repulsion. These terms are additional to those considered in 0, |^, which are 
based on the t:N interaction alone. Their derivation and interpretation within 
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the many-body scheme of this paper is given in [11|. In terms of an equivalent 
i^^o.ch parameter, using again the effective density of one obtains 

Sbo^^i, = -0.0022 ± 0.0002 m'^ (6) 

where the 10% error is an estimate due to second order effects in the modern 



chiral perturbation expansion |T^, which improves upon the Weinberg results. 

The present paper adds to these, the effects of short range nuclear cor- 
relations in conjunction with the off-shell extrapolation of the vrA^ scattering 
amplitude. Short range correlations are taken into account in the evaluation 
of the p-wave pion self-energy. Their effects are very important and lead to 
the famous Lorentz-Lorenz effect [H. The interpretation of this effect in terms 



of nuclear correlations modifying the ti + p exchange was given in [|I4|. In 
spite of its relevance in the modification of the p-wave self-energy, a similar 
work on the modification of the s-wave self-energy by the nuclear short range 
correlations has never been done and this is the purpose of this paper. We 



follow the ideas of [|14| to perform this work. 

In fig. 1 we show diagrammatically how a pion exchange in a rescattering 
term is modified due to the simultaneous exchange of other mesons (in the 
meson exchange picture) which generate the nuclear correlations. It is well- 



known [|T5| that the effect of this multiple exchange of mesons, simultaneously 
with the pion, can be approximated by means of a correlation function such 
that, if the diagram of fig. la can be represented by a potential V{r), the sum 
of la plus all those of the type lb can be represented by V{r), given by 

V{f) = V{f)9in (7) 

where 5'(r ) is an appropriate correlation function incorporating the repulsion 
at short distances. Hence 

(?(f) = l-/(f) with /(0) = 1, /(oo)=0 (8) 

Since in nuclear matter it is easier to apply Feynman rules in momentum 
space, we have instead 

-^^V{q-k)Q{k) (9) 

where 

V{k)= Jd^re'^''V{f) 

Q{k)= J d^re''''' g{f) = (27:^60) - B{k); B{k) = J d^re''''' f{f) 

(10) 

Hence 

V{q) = V{q) + K(g) ; K(g) = - / ^^(?"'- k)B{k) (11) 
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In this way we separate in V^(g ) one piece which corresponds to the ordinary 
one pion exchange and a second one induced by the correlations. 

If we particularize the diagrams of fig. 1 to the s-wave vrA^ interaction, 
we have the series for the pion propagator depicted in fig. 2. The terms 
connected by a one pion exchange line are summed up automatically by the 
Dyson equation, by writing the pion propagator which incorporates the proper 
pion self-energy. This self-energy corresponds to the irreducible diagrams (i.e. 
not connected by a pion line) in figs. 2b and 2d. There we see the lowest order 
s-wave pion self-energy (26) and a new piece generated by the correlations {2d) 
where the wavy line corresponds to V^(g ) of eq. (11). 

The interesting thing to observe is that, in the correlated potential, V 
enters as a function of q — k, with k an integration variable, and it is well- 
known that the isoscalar part of the s-wave nN amplitude (the one that enters 
in fig. 2d for isospin saturated nuclei) is small on-shell due to PCAC, but grows 
very fast as one moves to off-shell situations |]7|, |8|, p!6[ . 

In order to make a detailed evaluation of this piece let us recall that the 
s-wave ttN t-matrix (diagonal in spin) is written as [O 



[ 2Ai(g,gO ^ 



mt 



(12) 



with a, P, A pionic isospin indices. The lowest order pion self-energy corre- 
sponding to fig. 2b is written as 



n(g) 



2Ai(g,g) 
4tt p 



(13) 



and the second order term, induced by the correlations, corresponding to fig. 2d 
is given by 

1 2 



d^k 
(2^ 



2Ai(g,A;-g) 
An p 



{q - kf 



mt 



:B{k. 



(14) 



We observe that the isoscalar amplitude enters there as a function of the 
half off-shell variables. For pionic atoms we have q = [mj^jO) and q — k = 
(m^, —k), since we are assuming the exchange of mesons generating the cor- 
relations to be static (a fair approximation if they are of sufficiently short 
range) . 

In order to get a quantitative result for Hdq) with some estimation of the 
errors we have taken three off-shell extrapolations of the isoscalar amplitude 
and three different correlation functions. The three models used for the ex- 
trapolation are those of refs. 0, |[ and the recent one of the Jiilich group [|18 



The first one assumes for the isoscalar amplitude, of interest here, a short 
range piece together with a exchange and has been used in several many-body 



calculations P, |T9|, It gives, at g° = q 



J 



m^ 



Ai(g,g') = --(l + e)m^ 



m„ 



^1 + {q- q 



l\2 



(15) 
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with = 0.220 m-\ a,, = -0.233 m;^^ and = 550 MeV. 

The second one is based on a separable model for the vrA^ interaction and 
the use of dispersion relations. It can be parametrized in the region of interest 
to us as 1^ 



Ai(g, q') = Ai(0)(l + cip2 + + C3/) exp(-c/) (16) 

withp= \q-q'\, Ai(0) = 0.0075, Ci = 4.98 fm^, ca = 0.726 fm^ C3 = 0.203 fm^ 
and c = 0.462 fm^. 

The third model [|r^ uses (direct and crossed) nucleon and delta-isobar 
pole diagrams together with correlated 27r-exchange in the a (J=0) and p 
(J=l) channels. It describes quantitatively all relevant vrA^ scattering phase 



shifts over the whole elastic region. We take model (2) of [0 which, for 
reasons indicated in that paper, is preferred. The normalization of the T- 
matrix of |TB[ needed for our computation is the one of eq. (1) of [|T^ with 
the equivalence T*^^-* = 87rAi/m,r- In the region of interest to us this amplitude 
can be parametrized as 

Ai(g, q') = ih + hp + hp'') exp(-(p - po)V^') (17) 

with p = \q-q'\, h = 3.254 lO^^, = 0.1156 GeV^S h = 0.1588 GeV-^, 
Po = 0.204 GeV, and lj = 0.530 GeV. The quantity Ai is plotted in fig. 3 for 
the three models. 

On the other hand we have chosen three different correlation functions 
suited to calculations of nuclear matter: 

I) g{r) = 1 - joiq^r) 

II) g{r) = 1 - exp(-rV/?2) ^Ig) 
III) g{r) = (1 — exp(— r^/a^))^ + exp(— r^/c^) 

with = 780 MeV, R = 0.75 fm, a = 0.5 fm, b = 0.25 fm'^ and c = 1.28 fm. 

The results for are shown in Table I where we give the correction in 
terms of a coefficient SUbBq c- We can see that the contribution from this 
correlated piece is systematically repulsive, and the strength depends on the 
extrapolation and the correlations used. 

For a given off-shell extrapolation there is a certain dispersion of the results 
using different correlation functions. We use this to obtain an estimate of the 
uncertainties. Hence we take the average value and an error which allows us 
to reach the extremes. 

The results obtained with the Hamilton and Jiilich extrapolations are very 
similar, but they differ appreciably from those of the LMM model. Although 
we would favour the results of the Jiilich model, the large sensitivity to the 
off-shell extrapolation is an unpleasant feature. But we must recall here that 
the dispersion corrections of eqs. (5) are also dependent on the off-shell ex- 
trapolation. An interesting finding of the present work is the realization that 
there are cancellations between the two pieces which depend on the off-shell 
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extrapolation, the dispersion and the correlation corrections, and that when 
the terms are added, the sum is small compared to Bq + Abo and the depen- 
dence on the off-shell extrapolation is weak. This can be seen in Table II. The 
dispersive correction column has been completed with the results of the Jiilich 
model using the theoretical framework of |^ and the Jiilich extrapolation. 
For this purpose the extrapolation of A2 is also needed and we find for the 
region of interest to our problem 

A2(g,g') = A2(0)(l + ap')exp(-6p2) (19) 

with p = \q- q'\, A2(0) = 0.053, a = 9 GeV-^, b = 3.9 GeV-^. 

In the dispersive results of Table II we have included errors of 20% due 
to uncertainties other than the off-shell extrapolation p, p. The sum of the 
dispersive and correlation corrections appears in the third column of Table 
II and we observe that the results are: i) small, (of the order of 15% of the 
bo+Abo term), ii) rather similar for the different off-shell extrapolations and iii) 
compatible within the theoretical errors. We take an average of this sum and its 
errors as representative of this combined contribution and, after transforming 
it into a Sbo^d+c term, we show it below, together with the other contributions. 

If we add up all the contributions discussed in this paper, with their cor- 
responding uncertainties, to yield the theoretical equivalent value, b}^ we find: 

bo = -0.0130 ± 0.0033 lowest order 

A60 = -0.0140 ± 0.0008 m-i, second order Pauh 

corrected rescattering 

Sbo^d+c = +0.0039 ± 0.0020 m~^, dispersion + correlation ^'^^^ 

Sbo,ch = -0.0022 ± 0.0002 m-i, chiral 

= -0.0253 ± 0.0063 m^i, total 

The parameter Abo is evaluated from eq. (4) using again the effective density. 
The errors in bo and A60 are those coming from the experimental errors in ai 
and 03. 

It is interesting to compare these results to empirical values obtained from 
best fits to pionic atom data. Using the potential of Seki and Masutani 0, 
the one of Meirav et al. [^] and the one of Nieves et al. [22|, and using again 



the concept of effective density to write all of them in terms of an equivalent 
60 parameter, we find 

60= -0.0285 ± 0.0007 [9] 

bo= -0.0350 ± 0.0062 [21] (21) 

bo= -0.0323 ± 0.0007 [22] 



6 



where the errors are the statistical errors coming from the best fit (see section 



5 of ref. |23[ for a reanalysis with these potentials and the determination 
of the statistical errors). The results in eqs. (21) mean that there are also 
uncertainties in the empirical value of bo tied to different assumptions about 
the shape of the potential. It is also interesting to observe that although 
the theoretical result of eqs. (20) lies a little below the empirical results, it 
is compatible with all of them if uncertainties are considered. This becomes 
more obvious if a weighted average of the empirical results of eqs. (21) 



is done, which leads to b'^^ = —0.0304 ± 0.0005m~^, where again the error 
is only statistical. The dispersion of results in eqs. (21) indicates that there 
is an error of the order of 0.002m~^ tied to the assumptions in the shape of 
the potentials, such that the empirical value with its error becomes 6o™^ = 
—0.0304 ± 0.0025m~^. The central value of the theoretical results is 17% 
below the empirical one, but if uncertainties are considered the results are 
compatible. 

One last comment should be made about the role of the irNN form factor 
in the present calculations. We have used throughout a monopole form factor 
with a cut off mass A = 1300 MeV from the full Bonn A^A^ potential [^]. The 



form factor affects only the dispersion and chiral terms. There are indications 



that the cut off mass. A, could be smaller |^6[. In this case the strength of 
both terms is reduced by a similar percentage. Since the dispersive term is 
about three times larger than the chiral term, the net effect is to increase \&^\ 
and bring it even closer to the empirical value. 

The results of this paper have shown that correlation effects in the s-wave 
potential are sizeable and that the off-shell dependence of the nN amplitude 
must be considered in its evaluation. The correlation correction is then sensi- 
tive to the off-shell extrapolation used. However, when this correction is added 
to the dispersive correction, similarly sensitive to the off-shell behaviour, the 
sum of the two effects is small and rather insensitive to the extrapolation used. 

The present results should lead us to reflexion. In the past many theoretical 
efforts have been devoted to the search of the "missing repulsion" (see section 



4 of ref. [TT|] for a thorough and critical discussion of the different theoretical 
works). The results obtained here indicate that the "missing repulsion" , if any, 
cannot be large. They certainly do not exclude other sources of repulsion, like 
those tied to the renormalization of pions and nucleons in the inner structure 



of the vrA^ models [^, but puts severe constraints on their strength. 
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Table I : Correlation correction to the optical potential parameter KeBQlm'^] 
(eq. (1)) derived from the pion-nucleon amplitudes of Hamilton [0] , LMM [§ 
and the Jiilich group using three different correlation functions (eqs. (18)). 



SReBo^m-'] 


I 


II 


HI 


Hamilton 


-0.0127 


-0.0133 


-0.0141 


Jiilich 


-0.0051 


-0.0146 


-0.0094 


LMM 


-0.0001 


-0.0080 


-0.0012 



Table II: Dispersive and correlation corrections plus their sum in terms 
of the Rei?o parameter of eq. (1). The second column is the average of the 
correlation effects of Table I for the three correlation functions. 





SReBo^d 


5ReBo,c 


sum 


Hamilton [7] 


0.032 ±0.006 [5] 


-0.013 ±0.001 


0.019 ±0.007 


Jiilich [8] 


0.028 ±0.006 


-0.010 ±0.005 


0.018 ±0.011 


LMM [18] 


0.017 ±0.004 [6] 


-0.003 ±0.005 


0.014 ±0.009 
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FIGURE CAPTIONS: 

Fig. 1. (a) Second order term in nA scattering mediated by pion exchange 
alone. (6) Schematic representation of simultaneous exchange of other mesons 
together with pion exchange, which are taken into account in terms of a cor- 
relation function. 

Fig. 2. Diagrammatic representation of the pion propagator led by the pion 
s-wave self-energy. The piece (d) is mediated by the correlated potential of 
eq. (11). 

Fig. 3. The value of Xi{0,k) defined in eqs. (15), (16) and (17) as a function 
of the momentum transfer in units of pion mass, for pionic atom kinematics. 



(a) Hamilton model from ref. [0]. (b) Jiilich group model (2) from ref. [|T8 
(c) LMM model from ref. HI. 
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Fig. 2 




Fig. 3 



